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Abstract 

The antiferromagnetic Heisenberg model on a pyrochlore lattice under external magnetic field is studied by 
classical Monte Carlo simulation. The model includes bilinear and biquadratic interactions; the latter effectively 
describes the coupling to lattice distortions. The magnetization process shows a half-magnetization plateau at low 
temperatures, accompanied with strong suppression of the magnetic susceptibility. Temperature dependence of 
the plateau behavior is clarified. Finite-temperature phase diagram under the magnetic field is determined. The 
results are compared with recent experimental results in chromium spinel oxides. 
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1. Introduction 

Geometrically frustrated magnetism is one of tlie 
long-standing issues in condensed-matter physics. 
Geometrical frustration leads to nearly-degenerate 
ground-state manifolds of a large number of differ- 
ent spin configurations. The degeneracy may yield 
nontrivial phenomena, such as spin liquid states 
and glassy states [1,2]. Because of the small energy 
scale of the nearly-degenerate manifolds, a tiny 
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perturbation can be relevant to lift the degeneracy 
in these systems. For example, quantum or ther- 
mal fluctuations can give rise to non-trivial phase 
transitions in applied magnetic field, leading to a 
plateau or jump in the magnetization. 

Pyrochlore lattice is a typical example of the 
geometrically-frustrated structures; it consists of 
a three-dimensional network of corner-sharing 
tetrahedra as shown in Fig. 1. Pyrochlore lattice 
is highly frustrated: Classical antiferromagnetic 
Heisenberg models with only nearest-neighbor in- 
teractions do not show any long-range ordering 
down to the lowest temperature T = [3,4]. The 
situation is similar in the S — 1/2 quantum spin 
model as well [5,6]. Effects of several perturba- 
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Fig. 1. Cubic unit cell of the pyrochlore lattice. Near- 
est—neighbor exchanges J and third— neighbor exchanges 
J3 are shown, which are incorporated in our model for the 
MC calculations. See Sec. 2 for the details. 

tions on the highly-degenerate states have been 
discussed, such as quantum fluctuations [6,7], 
further-neighbor interactions [3,8], spin-lattice 
coupling [9,10] and Dzyaloshinsky-Moriya inter- 
action [11]. 

Pyrochlore magents are found in many real com- 
pounds. Most typically, so-called B spinel oxides 
ABiO^ have the pyrochlore structure of the mag- 
netic B cations. Among the B spinels, chromium 
spinel oxides ACr2 04 with nonmagnetic A cations 
have attracted much interests recently [12]. In 
these compounds, each Cr'^^ cation has three d 
electrons in three-fold t2g levels, which constitute 
the S — 3/2 localized spin by the Hund's-rule 
coupling. Hence, t^g levels are half-filled and there 
is no orbital degree of freedom in contrast to the 
compounds with B=N [13,14] or Ti [15]. There- 
fore, the Cr spinels provide simple S = ij^ spin 
systems where detailed comparisons between ex- 
perimental and theoretical studies can be made 
rather straightforwardly. 

Recently, the chromium spinel oxides have been 
studied by applying external magnetic field [16]. 
Surprisingly, it was revealed that some of them ex- 
hibit a half-magnetization plateau in a very wide 
range of the magnetic field. For instance, the mag- 
netization process of HgCr204 shows a plateau 
from 10 Tesla to 27 Tesla. 

The spin-lattice coupling as a possible mecha- 



nism of the plateau formation in the Cr spinels has 
been proposed and studied theoretically [17]. The 
spin-lattice coupling leads to effective biquadratic 
interaction between spins. The T = study of the 
classical Heisenberg model on the pyrochlore lat- 
tice has revealed that for finite biquadratic inter- 
actions, the model exhibits a coUinear state at in- 
termediate magnetic-field regime, in which each 
tetrahedron has 3-up and 1-down spin configura- 
tion. This coUinear state gives rise to a stable half- 
magnetization plateau. The magnetization process 
is favorably compared with the experimental re- 
sults. 

Main purpose of the present study is to investi- 
gate the plateau problem in the pyrochlore Heisen- 
berg model at finite temperatures in order to make 
further comparisons with the experimental results. 
We will employ Monte Carlo (MC) calculations to 
clarify effects of thermal fiuctuations, which are ne- 
glected in the previous T = study and may play a 
crucial role in the frustrated systems. We will show 
the temperature dependences of the magnetization 
curve and the magnetic susceptibility under the 
external magnetic field, and summarize the finite- 
temperature phase diagram. We will compare the 
results with the experimental data. 

This paper is organized as follows. In Sec. 2, 
we introduce the bilinear-biquadratic Heisenberg 
model. MC method is also described. In Sec. 3, 
we show numerical results in comparison with the 
previous T = results and the experimental data. 
Section 4 is devoted to summary and concluding 
remarks. 

2. Model and Method 

Our starting point is a classical Heisenberg 
model with the spin-lattice coupling of spin- 
Peierls type on the pyrochlore lattice under the 
external magnetic field [17]. The Hamiltonian is 
given in the form 

-~Y.h-s., (1) 



where Jij is the exchange constant for the spins at 
the sites i and j, S is a classical vector with the 
modulus |S| = 1, a is the spin-lattice coupling 
constant, pij is the change of the bond length be- 
tween «th and jth sites, Kij is the elastic constant 
for the ij bond, and h is the external magnetic 
field. This model has been studied at zero magnetic 
field to describe the low-temperature properties of 
Cr spinel compounds [10]. 

In Eq. (1), the lattice distortions pij are 
quadratic, leading to Gaussian integrals in the 
partition function, and therefore they can be in- 
tegrated out. In general case, the integration re- 
sults in a complicated effective spin-Hamiltonian 
with long-range interactions between spin-bonds. 
However, assuming four-sublattice long-range 
ordering, and nearest-neighbor exchange J and 
elastic coupling K only [17], the effective Hamil- 
tonian for the ordered moments appearing in the 
resultant partition function is equivalent to 

H = Y, J [Si ■ Sj - b{Si ■ Sjf} -J2h-Si, (2) 

{ij) » 

where the summation is taken over the nearest- 
neighbor sites. The second term in Eq. (2) is the 
biquadratic interaction, with the coupling con- 
stant b — Ja^/K. Hence, this biquadratic inter- 
action mimics the spin-lattice coupling. It favors 
a collinear spin configuration as the spin-lattice 
coupling does since Jb is always positive. 

In the present study, for simplicity and as a first 
step, we neglect all the complications which arise 
from the lattice distortions apart from the effective 
bilinear- biquadratic part, and consider the model 
(2) extended with the third-neighbor exchange J3 
(as shown in Fig. 1). In particular, we focus on the 
ferromagnetic J3 < in the following calculations 
since it stabilizes the simple four-sublattice long- 
range ordering with the wave vector q = [3]. 
Hereafter, we will set J = 1 as an energy unit 
and use the convention of the Boltzmann constant 
kB = 1. 

We will investigate thermodynamic properties of 
the model (2) by Monte Carlo (MC) calculations 
[18] . We employ the Metropolis algorithm with lo- 
cal spin updates to sample spin configurations. We 
typically perform 10^ MC samplings for measure- 



ments after 2 x 10^ steps for thermalization. The 
measurements are performed in every iVint-times 
MC update, and we typically take A^int = 2. Re- 
sults are divided into five bins to estimate statisti- 
cal errors by variance of average values in the bins. 
We have checked the convergence and ergodicity of 
the results by comparing those for different initial 
spin configurations. The system sizes in the present 
work are up to i = 8, where L is the linear di- 
mension of the system measured in the cubic units 
shown in Fig. 1, i.e., the total number of spins N 
is given by x 16. 

3. Results and Discussion 

First, we show the magnetic properties of the 
model. Figure 2 shows the magnetic field depen- 
dences of the magnetization and the uniform mag- 
netic susceptibility at different temperatures for 
J3 = —0.05 and b = 0.1. The magnetization is de- 
fined as 

m=(Mtot) = (|^5i|/7v), (3) 

i 

and the susceptibility is measured by 

X = N{{Ml,)-{M,,,f)/T, (4) 

where the brackets denote the thermal average. We 
show the data for L = 8 since there is no significant 
system size dependence in the range of parameters 
investigated. 

As shown in Fig. 2, at low temperatures, the 
magnetization shows a clear plateau feature at 
around the half magnetization m = 1/2. At the 
same time, the susceptibility exhibits a sharp dip 
in the corresponding range of the magnetic field. 
The results indicate that the plateau is stable and 
survives up to T ~ 0.3J for the present parame- 
ters. 

Next, we examine the phase diagram under the 
magnetic field. Following Rcf. [17], we measured 
the order parameters in the form 

><. = J2J2<c/NteU., (5) 

tetra a 

where the first summation is taken over all 
the Ntetia. = -^/4 independent tetrahedra (up- 
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Fig. 2. Monte Carlo data at finite temperatures of the 
magnetic-field dependences of (a) the magnetization and 
(b) the uniform susceptibiUty for J3 = —0.05 and h = 0.1. 
The Hnes are guides for the eyes. The bold gray curves 
show the results at T = [17]. 

pointing ones in the [111] direction). The index 
v denotes the irreducible representations (irreps) 
of the tetrahedral symmetry group T^, i.e., v = 
{Ai , E, T2}, and the summation on a is taken over 
members of a given irreps. The irreps A^.q, are 
given by p^^a in Eq. (4) in Ref. [17] with replacing 
Pij by Si ■ Sj on r.h.s of the equation. 

Figure 3 shows the MC results for the order pa- 
rameters with the symmetry E and T2. At low 
temperatures, Xe becomes finite in the low-field 
regime h ^ 0. This corresponds to a coplanar long- 
range ordering with tetragonal symmetry. On the 
other hand, becomes finite in a middle range of 
the magnetic field h ^ 4, indicating a 3-up 1-down 



Fig. 3. Monte Carlo data of the order parameters with 
(a) E symmetry and (b) T2 symmetry for J3 = —0.05 and 
b = 0.1. The lines are guides for the eyes. The bold gray 
curves show the results at T = [17]. 

type long-range ordering with trigonal symmetry. 
Typical spin configurations are shown in Fig. 4. 

The order parameter alone, however, does 
not differentiate between the coUinear phase which 
exhibits the magnetization plateau in Fig. 2 and 
the non-coUinear but coplanar spin-canting phase 
as shown in the phase diagram at T = [17]. 
To identify these phases, we have examined the 
coUinearity and the coplanarity of the spins by 
calculating the so-called nematic order parame- 
ters, and successfully determined the phase bound- 
ary between the two phases at finite temperatures 
[19]. We note that the obtained coUinear-coplanar 
phase boundary corresponds to the shoulder-like 
feature in the susceptibility at /i ^ 4 — 5 in Fig. 2 
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Fig. 4. Monte Carlo phase diagram for J3 = —0.05 and 
b = 0.1 Symmetry and typical spin configuration in the 
tetrahedron are shown for each phase. The shaded area 
denotes the half-magnetization plateau phase. Filled sym- 
bols are the T = results [17]. The curves for the phase 
boundaries are guides for the eyes. 

Figure 4 summarizes the MC phase diagram for 
J3 = —0.05 and b = 0.1, determined from the data 
for the system sizes L — 4,6 and 8. All the phase 
transitions are of first order, and the phase bound- 
aries smoothly extrapolate the zero-temperature 
theoretical values. In addition to the already dis- 
cussed three ordered phases — the E phase, the 
collinear T2 phase, the spin-canted T2 phase — at 
high temperatures there is the disordered param- 
agnetic phase. In the figure we also show the typi- 
cal four-sublattice spin configurations on the tetra- 
hedron for the three ordered phases. The plateau 
phase is robust against finite-temperature fluctua- 
tions, and survives up to even slightly higher tem- 
perature than the critical temperature at zero mag- 
netic field. 

The present MC results can be compared favor- 
ably with the recent experimental results in Cr 
spinel compounds such as CdCr2 04 and IIgCr204. 
In particular, the magnetization process in the 
HgCr2 04 has been measured up to the saturation 
field and the entire phase diagram has been deter- 
mined [16]. Our MC results in Figs. 2 and 4 give a 



qualitative description of the experimental results. 
More quantitative comparison including a fine 
tuning of the model parameters will be reported 
elsewhere [19]. 

4. Summary and Concluding Remarks 

To summarize, we have investigated the mag- 
netic properties at finite temperatures of the anti- 
ferromagnetic Heisenberg model on the pyrochlore 
lattice by Monte Carlo calculations. In addition 
to the ordinary bilinear exchange interaction, our 
model includes the biquadratic interaction as an 
effective term describing the coupling of the spins 
to lattice distortions. By applying magnetic field, 
we have shown that the model exhibits a stable 
half-magnetization plateau in a considerable range 
of temperatures. The results show good agreement 
with recent experimental results in Cr spinel ox- 
ides. 

Our study indicates that the spin-lattice cou- 
pling plays a dominant role in stabilizing the 
half-magnetization plateau in Cr spinel oxides. 
However, we note that other additional elements 
are indispensable to understand the properties of 
real compounds in more detail. For instance, at 
zero magnetic field, recent neutron experiments 
have found that the spin ordering in the low- 
temperature phase is complicated, compound- 
dependent, and incommensurate [20]. To account 
for such effects, further studies of more realistic 
Hamiltonians are needed, which would include, 
e.g., longer-range exchange interactions. 

Our further calculations for the model (2) sug- 
gest an interesting phenomenon when we weaken 
the strength of the biquadratic interaction b. For 
small values of b, the width (along the h axis) of 
the plateau phase grows with the temperature. 
This clearly implies that thermal fluctuations 
also stabilize the plateau phase. This tendency 
becomes conspicuous as b decreases, opening a 
route toward the so-called order-from-disorder 
phenomenon [21]. This interesting issue will be 
discussed elsewhere [22]. 

Finally, our model is defined for classical spins 
although the real compounds have S = 3/2 quan- 
tum spins. We have also studied the quantum spin 



problem by spin wave analysis, and found that, as 
expected, the plateau state has a finite spin gap. 
The details of the spin excitation spectrum will 
also be discussed elsewhere [23] . 
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